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Abstract 



Using the CLEO II detector at the Cornell Electron Storage Ring we have 
measured the H|? decay asymmetry parameter in the decay H[! —>■ H~7r^. We 
find a^oa- = 0.26 ± 0.18(stat)+o;o4(syst), using the world average value of 
OH = -0.456 ± 0.014 we obtain a^o = -0.56 ± 0.39(stat)+°;o°(syst). The 
physically allowed range of a decay asymmetry parameter is— l<a<+l. 
Our result prefers a negative value: a=o is < 0.1 at the 90% CL. The central 
value occupies the middle of the theoretically expected range but is not yet 
precise enough to choose between models. 

PACS numbers 14.20.Lq,14.20.Jn,14.65.Dw,11.30.Er 
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The weak interactions underlying hadronic charm quark weak decay are straightforward 
to describe theoretically, but complications arise because the quarks are bound inside hadrons 
by the strong force. These interactions, which are described by the theory of quantum 
chromodynamics (QCD), are very difficult to predict using perturbative methods because the 
strong coupling is large at the typical energies of charm decays. Compared to charm mesons, 
charm baryons offer new information for two reasons: first, non-factorizable contributions 
to the decay amplitude are important; W-exchange diagrams can contribute without the 
helicity suppression that decreases their contribution to pseudoscalar meson decays, and 
internal W emission is expected to be significant. The relative importance of these effects is 
believed to be responsible for both the observed lifetime hierarchy of the charm baryons and 
the relatively shorter lifetime of charm baryons compared to charm mesons. Second, parity 
violation in charm baryon decays is readily observable because the decay of the daughter 
hyperon also violates parity. In consequence a variety of models ^ are used to predict 
both the decay rate and degree of parity violation in charm baryon decays. To constrain the 
models it is important to provide as much experimental information as possible. 

Parity violation occurs in hyperon and charm baryon |^ ^^0 decays due to the 
existence of two orbital angular momentum decay amplitudes of opposite parity. The ex- 
perimental observable is an asymmetry in the angular decay distribution due to interference 
between the two amplitudes. In the decays — > Att^ and S~ — > Avr" both followed by 
A pTT~ , the A is produced with a polarization equal to 

^ + A-Pb)A-/?b(A X Pb) -7bA X (A X Pb) 

"a — ' ( J- ) 

(l + a^A-PB) 

where Pb is the parent baryon polarization, ob, Pb and 75 are the parent baryon asymmetry 
parameters and A is a unit vector along the A momentum in the parent baryon frame [0. 
If the parent baryon polarization is unobserved, or if the parent baryon is not polarized. 
Equation |l] reduces to Pa = A. The angular distribution of a proton from the decay of 
a A is therefore 

dN 

— oc 1 + ttsaACOsfc'A (2) 

a cos b'A 

where ^a is the angle between the proton momentum vector in the A rest frame and A, 
and ttA is the A decay asymmetry, which measures the degree of parity violation in A decay. 
Experimentally, a fit to equation ^ determines a product of decay asymmetries. For example, 
in Att~, as = since is well measured (q;a = 0.642 ± 0.013 ), a-~- can be 

extracted. 

Equation ^ has also been used to extract the decay asymmetry of the A^ in several decay 
modes. For example in A+ Att"*", = In the massless fermion limit the chirality 

of the weak interaction corresponds to a = —1, however, the interplay between the strong 
and weak interactions in the weak decay alters the value of a. The decay asymmetry in 
A+ Atv^ is relatively well measured, the world average value is = —0.98 ± 0.19 [Q, 
which is consistent with the naive V-A expectation. This result had been predicted by Heavy 
Quark Effective Theory (HQET) in conjunction with the factorization hypothesis [Q. Most 
pole and quark models [0] , also predict large negative values of q:^+ . 
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In this letter we measure the decay asymmetry parameter in the decay S^7r+ 
for the first time. For — > S^tt^ theoretical predictions cover a larger range: —1 < a^o < 
—0.38 The decay of S° is a three step process where parity violation occurs at each 
decay stage: S° — > H^7r+, — > Ayr", A — >• p7r~. The differential rate is therefore given 
by §: 

dN 

— ; — oc 1 + a=oa= cos 6^= + a^a^^ cos 

a cos tt cos t^A 

+ aA«HO COS 6's cos ^A (3) 

where in addition to the terms previously defined: 6^ is the angle between the A momentum 
vector in the S~ rest frame and the S~ momentum vector in the S° rest frame as shown in 
Fig. |l[ Equation ^ reduces to the familiar form of Equation || when integrated over cos ^a- 
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FIG. 1. Definition of the angles 0= and ^a in S[; ^ H vr^. The decay of each particle is drawn 
in its own rest frame. 

The data sample in this study was collected with the CLEO II detector [^] at the Cornell 
Electron Storage Ring (CESR). The integrated luminosity consists of 4.83 fb^^ taken at and 
just below the T(4iS) resonance, corresponding to approximately 5 million e~^e^ cc events. 

We search for the decay S° S^tt^ in e+e^ cc events by reconstructing the decay 
chain S° S~7r^, Att^, A Ptt^Q. 

The A is reconstructed by requiring two oppositely charged tracks to originate from a 
common vertex. The positive track is required to be consistent with a proton hypothesis^. 
The momentum of the A candidate is calculated by extrapolating the charged track momenta 
to the secondary vertex. The invariant mass of A candidates is required to be within three 
standard deviations (3cr = 6.0 MeV/c^) of the known A mass. Track combinations which 



Unless otherwise noted, throughout this letter charge conjugation is implied 

^Hadronic particles are identified by requiring specific ionization energy loss measurements 
(dE/dx), combined with time-of- flight (TOF) information when available. The two measurements 
are combined into a joint probability for the particle to be a pion, a kaon or a proton. A charged 
track is defined to be consistent with a particle hypothesis if its probability is greater than 0.003. 
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satisfy interpretation as tt+tt" are rejected. Combinatoric and B decay backgrounds 

are reduced by requiring the momentum of A candidates be greater than 800 MeV/c. 

The S~ is reconstructed in the Avr" decay mode. S~ candidates are formed by combining 
each A candidate with a negatively charged track consistent with a pion hypothesis. The 

candidate vertex is formed from the intersection of the A momentum vector and the 
negatively charged track. To obtain the S~ momentum, and An invariant mass, the momen- 
tum of the charged track is recalculated at the new vertex. The invariant Avr" mass is fit 
to a double Gaussian signal shape with parameters fixed from a GEANT based Monte 
Carlo (MC) simulation of the detector, and a first order Chebyshev polynomial to describe 
the combinatorial background. We find 8434 ± 109 events consistent with Att~ with 

an invariant mass of 1322.0 ± 0.03 MeV/c^ and width = 2.5 MeV/c^, where is the 
weighted average width of the double Gaussian. The mean and width are in agreement with 
the MC simulation. 

To reconstruct H° candidates all S~ candidates are combined with positively charged 
tracks consistent with a pion hypothesis. Since charm fragmentation is a relatively hard 
process, the S° spectrum from e~^e~ cc is also fairly stiff. We therefore use a scaled 
momentum cut, Xp > 0.5 Q to reduce combinatoric background. This cut also eliminates 
charmed baryons from decays of B mesons. Fig. § shows the invariant mass distribution 
of S~7r"'" combinations. The S~7r"'" invariant mass distribution is fit to a double Gaussian 
to describe the signal and a first order Chebyshev polynomial to describe the combinatorial 
background, where the parameters of the double Gaussian are fixed by MC simulation. 
The shaded histogram is wrong sign (WS) random S~7r~ combinations. There are less WS 
than right sign (RS) random combinations due to charge conservation. The excess of right 
sign events over wrong sign events below M{'E~n~^) = 2.33 GeV/c^ is due to feedthrough 
from S^n(7r), n > 1, where one or more pions are missing. The kinematic limit 

for feedthrough is M(S°) - M{n) = 2.3304 GeV/c^. To simplify the fit to M(S-7r+) the 
feedthrough region is excluded. We find 138 ± 14 ^ H^7r+ candidates with a mass of 
2470.6 ± 1.0 MeV/c^ and width aav = 10 MeV/c^. The mean and width are in agreement 
with MC simulation. We require that the invariant mass of H° candidates be within 30 
MeV/c^ (3 aav) of the known S° mass. 

MC simulation shows that the difference between generated and reconstructed values 
of cos^H and cos6'a has both Gaussian and symmetric non-Gaussian components. The 
resolution in cos 6'= a, cr(cos 6'=, A), defined to be the average of the rms variance of the 
Gaussian and of the non-Gaussian component, weighted by the relative normalizations of 
the two components Q is: cr(cos6's) = 0.014 and cr(cos6'A) = 0.030. 

The decay asymmetry parameter in — * S^tt^ is measured in a two dimensional un- 
binned maximum loglikelihood fit to the two-fold decay angular distribution of Equation ^ 
in a manner similar to This technique enables a multi-dimensional likelihood fit to be 



■^For cos^H and cos6'a, the non-Gaussian component comprises 60% and 32% of the distribution, 
respectively 
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FIG. 2. Invariant H 
vr^ combinations. 



2.2 2.3 2.4 2.5 2.6 2.7 

7r+ mass distribution. The shaded histogram is wrong sign (WS) random 



performed to variables modified by experimental acceptance and resolution. The probability 
function of the signal, F^, is determined by generating one high statistics MC sample of 
S° E~7i~^, E~ Att", a p7T~ with a known value of a^^o and the world average values 
of cts and a\. The generated events are processed through the detector simulation, off-line 
analysis programs, and selection criteria. Using the generated angles, accepted MC events 
are weighted by the ratio of the decay distribution for a trial value of a^o to that of the gen- 
erated distribution. By such weighting, a likelihood may be evaluated for each data event 
for trial values of a^o, and a fit performed. The probability for each event is determined by 
sampling Fg using a search area centered on each data point. The size of the area is chosen 
such that the systematic effect from a finite search area is small and the required number of 
MC events is not prohibitively high. 

Background is incorporated into the fitting technique by constructing the log-likelihood 
function: 

N 

lnC = Y.lniPsTs + PBTB) (4) 

i=l 

where is the number of events in the signal region and Ps and Pb are the probabilities that 
events in this region are signal and background respectively. The probability distribution of 
background in the signal region, F^, is determined from (H^tt^) mass sidebands above and 
below the signal region. The sidebands are: 2.5004 < M(E:-7r+)<2.6904 GeV/c^ and 2.3304 
< M(S-7r+)<2.4404 GeV/c^. 

CP conservation requires a^o = — a=o. As the number of signal events is small the analysis 
is insensitive to the presence of CP violation, therefore CP is assumed to be conserved. Since 
Equation ^ depends on the products a^^oaE., a^^oa^ and q:s«a which have the same sign for 
particle and anti-particle states, particle and anti-particle distributions are combined. 

The validity of the analysis procedure is determined by constructing artificial data sets 
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consisting of MC signal events, generated with known asymmetries, and background events 
taken from data. The generated asymmetry was varied over the range: — 1 < < 0, in all 
cases the fit returns unbiased values of a^o. 

Applying the fit procedure to the data we find asocts = 0.26 ±0.18, a^a\ = — 0.57±0.21 
and a^oc^A = —0.26 ±0.37. The product a^oa\ is second order in decay angles, consequently 
the fit is less sensitive to this quantity than to a^oa^ and a^a\. We provide a value of a^oa\ 
for completeness only. The cos 6^ and cos 6\ distributions for data and projections of the fit 
are shown in Fig. 0. The value of a^a\ found in the fit is in reasonable agreement with a 




FIG. 3. Upper plot: the cosine of the angle, 6^, between the A momentum vector in the H~ 
rest frame and the H~ momentum vector in the rest frame, in the decay sequence Hj? — > H^vr'^ 
H~ — > Avr", A — > pTr~ for sideband subtracted data (points with error bars) and projection of 
the fit (solid histogram). Lower plot: the cosine of the angle, 6\, between the proton momentum 
vector in the A rest frame and the A momentum vector in the H~ rest frame in the decay sequence 
H[! H^vr^ Avr^, A — > pvr" for sideband subtracted data (points with error bars) and 

projection of the fit (solid histogram). 

CLEO measurement of a^a^ = — 0.291±0.021(stat)^Q'pj^5(syst) obtained with approximately 
8,000 S~ — > Att" events and with the current world average value a;s«A = —0.293 ±0.007 
0. Constraining a^a\ to the world average value does not reduce the statistical error on 

c 

We have considered the following sources of systematic error and give our estimate of 
the percentage error on 0=0 in parentheses. The statistical error in the MC sample is 
estimated by varying the size of the MC sample used in the fit (±2%). The error associ- 
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ated with the uncertainty in the fragmentation function [|T^] ^ is estimated by varying 
this function (+1%, —14%). To determine the effect of incomplete knowledge of the back- 
ground shape and the effect of statistical fluctuations in the sideband sample used to model 
the background in the signal region, we vary the size of the lower and upper sidebands 
used in the loglikelihood fit. We also repeat the fit using wrong sign S^tt" events in both 
the signal mass region and sideband regions to model the background shape (+16%, — 1%). 
The error associated with MC modeling of slow pions from S and A decay is obtained by 
varying the reconstruction efficiency according to our understanding of the CLEO II detec- 
tor (+9%, —6%). Using a large sample of S° — > H~7r+ MC generated with q;=o = —0.5 
and the world average values of and and including a randomly generated back- 
ground based on the shape of the data sideband, we measure the effect of varying the 
size of the area element used to determine r5(cos6's, cos6'a; ol-o, c^a) and the background 
shape rB(cos6'H, cos^a) in the loglikelihood fit (+1%,— 2%). The fit method is checked 
by integrating over cos^a (cos^s) and performing a one dimensional binned fit to cos^s 
(cos^a)- The results are consistent with the maximum likelihood fit. Possible background 
from S°(2450) A+tt^j^^, A+ — * ^'^tast determined to be negligible by MC simulation. 
This measurement is insensitive to production polarization, Psg? and no systematic error 
has been included from this source [l^. Adding all sources of systematic error in quadra- 



ture (+19%,— 15%), we find a=oas = 0.26 ± 0.18(stat)^gQ4(syst). A similar study for the 



0.04V 

systematic uncertainty on the measurement of asctA in S° S^tt^ yields (—8%, +4%) and 



consequently asctA = — 0.57 ± 0.21(stat)^Qg2(syst). 

From the measurement of o;soa= and the PDG evaluation: a= = —0.456 ± 0.014 [Q, 
we obtain a=o = —0.56 ± 0.39(stat)^Q Q9(syst). The physically allowed range of a decay 
asymmetry parameter is — 1 < a < +1. Our result prefers a negative value: aco is < 0.1 
at the 90% CL. The central value is in the middle of the theoretically expected range but is 
not yet precise enough to choose between models. 

We note that has, so far, always been determined as a product of decay asymme- 
tries using Equation ^. In principle, given sufficient statistics, the loglikelihood method 
in conjunction with the two-fold joint angular distribution of Equation |^ allows the direct 
measurement of all three asymmetry parameters: a=o, 0;=, and ctA- 

In conclusion, from a sample of 138 ± 14 S° ^ H^tt^ decays we have measured 0=00= = 
0.26 ± 0.18(stat)^Q Q4(syst), from which we obtain = —0.56 ± 0.39(stat)^Q gg(syst). To 
the best of our knowledge this is the first measurement of a charm strange baryon decay 
asymmetry parameter. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, Research Corporation, the Natural Sciences and 
Engineering Research Council of Canada, the A. P. Sloan Foundation, the Swiss National 
Science Foundation, and the Alexander von Humboldt Stiftung. 



^Based on a measurement of the fragmentation function at CLEO, we use a Peterson fragmen- 
tation function |11] with = 0.15. 
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